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Abstract. We use time-frequency methods for the study of Fourier Integral 
operators (FIOs). In this paper we shall show that Gabor frames provide very 
efficient representations for a large class of FIOs. Indeed, similarly to the case 
of shearlets and curvelets frames I27| , the matrix representation of a Fourier 
Integral Operator with respect to a Gabor frame is well-organized. This is used 
as a powerful tool to study the boundedness of FIOs on modulation spaces. As 
special cases, we recapture boundedness results on modulation spaces for pseudo- 
differential operators with symbols in M 00 ' 1 [25], for some unimodular Fourier 
multipliers [2] and metaplectic operators [TQ1 H3] • 



Fourier Integral Operators (FIOs) are a mathematical tool to study variety of 
problems arising in partial differential equations. Originally introduced by Lax 
[55] for the construction of parametrices in the Cauchy problem for hyperbolic 
equations, they have been widely employed to represent solutions to Cauchy prob- 
lems, in the framework of both pure and applied mathematics (see, e.g., the papers 
[21 El [121 H31 [271 EQ], the books [3H [351 ES] and references therein). In particular, 
they were employed by Helffer and Robert [281 129] to study the spectral prop- 
erty of a class of globally elliptic operators, generalizing the harmonic oscillator of 
the Quantum Mechanics. The Fourier Integral operators we work with, possess a 
phase function similar to those of [23, EH] • A simple example is the resolvent of the 
Cauchy problem for the Schrodinger equation with a quadratic Hamiltonian. 

For a given function / on M. d the Fourier Integral Operator (FIO) T with symbol 
a and phase $ on ~R 2d can be formally defined by 
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1. Introduction 



Tf(x)= [ e^^M^Hv)^- 
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and the nondegeneracy condition 

(2) | det 8% >v $(x, g)\>5> 0, (x, rj) e K M 
The symbol a on M. 2d satisfies 

|<9>(z)| < C a , \a\ < N, a.e. z G R 2d , 

for a fixed > (in the sequel we shall work also with rougher symbols). 

The first goal of this paper is to rephrase the operator T in terms of time- 
frequency analysis (see Grochenig [23] and the next Section 2 for a review of the 
time-frequency methods.) Denoting T x f(t) = fit — x), M v f(t) = e 2mr,t f(t), for 
a, /3 > 0, g G L 2 (R d ), the set of time-frequency shifts Q(g,a,/3) = {g m) n '■ = 
M n T m g} with (m,n) e aZ d x(3Z d , is a Gabor frame if there exist positive constants 
A, B > 0, such that 

(3) A\\f\\ L ,<J2\(f,T m M n g)\ 2 <B\\f\\ L2 , V/GL 2 (M rf ). 

m,n 

In Section 3 we show that the matrix representation of a FIO T with respect to a 
Gabor frame with g e iS(M d ) is well-organized (similarly to frames of curvelets and 
shearlets [H1EZ1), provided that the symbol a satisfies the decay estimate for every 
N > (see Theorem 13. 3ft : 

Theorem 1. For each N > 0, there exists a constant > such that 

(4) \{Tg m! n,g m >,n>)\ < C N (x(m,n) - (m ,n'))~ 2N , 
where \ is the canonical transformation generated by 

In the special case of pseudodifferential operators such an almost diagonalization 
was already obtained in [2H [31] . Indeed, notice that pseudodifferential operators 
correspond to the phase $(x, g) = xg and canonical transformation x{Ui v) = {Vi v)- 

As a rilevant byproduct of the results of Section 3, we study the boundedness 
properties of the operator T on the so-called modulation spaces (Section 4 and 5). 
To define them, we fix a non-zero Schwartz function g and consider the short-time 
Fourier Transform V g f of a function / on M d with respect to g 



Vgfix, g) = (f, M v T x g) = / f(t) g(t - x) e~ 2 ^ dt 



which provides a time-frequency representation of /. The (unweighted) modulation 
space M p ' q is the closure of the Schwartz class with respect to the norm 

MP , q = \\V 9 f\\ LP , q = ( [ ([ \V g f(x, V )\ p dxY P d V 
\Jm d \JR d J 
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(with appropriate modifications when p = oo or q = oo). In particular, when 
p = q we simply write M p ' p = M p , see Subsection 12.21 for exhaustive definitions 
and properties. 

These spaces were introduced by Feichtinger [17] and have become canonical for 
both time-frequency and phase-space analysis [18J , most recent employment being 
the study of PDEs [21 El M> GSJ HQ] . 

If g G M 1 , and the Gabor frame {T m M n g; (m, n) G aZ d x f3Z d } is a tight frame, 
namely holds with A = B, then it extends to a Banach frame for the modulation 
spaces M p,q (M. 2d ), with the norm equivalence 

||/||jwp.4 x \\{f,T m M n g) mtn \\p>,q. 

Then, boundedness of the FIO T on M p,q is equivalent to that of the infinite matrix 
(Tg mt n, 9m', n') on the spaces of sequences l p,q . 

Whence, the estimates fll]) readily give (see Theorem 14.11 for a more general 
version) : 

Theorem 2. For N > d, 1 < p < oo, the Fourier integral operator T, with 
symbol a and phase $ as above, extends to a continuous operator on M p . 
(In the case p = oo, the space M°° is replaced by the closure of the Schwartz 
function with respect to || • \\m°° )■ 

The continuity property of a FIO T on M p,q , with p ^ q, fails in general. Indeed, 
an example is provided by the operator Tf(x) = e m \ x \ f(x), corresponding to 
rj) = XT] + |x| 2 /2, cr = l, which is bounded on M p ' q if and only if p = q (see 
Proposition 17.11) . 

Hence, we introduce a new condition on the phase $, namely that the map 
x i — > V x $(x,?7) has a range of finite diameter, uniformly with respect to rj, that 
allows us to get the boundedness on M p,q (Theorem 15 . 2j) : 

Theorem 3. For N > d, 1 < p,q < oo, the Fourier integral operator T , with 
symbol a and phase $ as above, and such that 

sup |V x .$(x,r?) - V x $(x',ri)\ < oo, 

x,x' ,r) 

extends to a continuous operator on M p,q . 

(In the case p = oo or q = oo, the space M p ' q is replaced by the closure of the 
Schwartz function with respect to \\ ■ \\mp<i)- 

As a particular case, we recapture recent boundedness results of unimodular 
Fourier multipliers [2] (see Example 15.31) . 

With respect to Theorem 2, here the proof is more delicate and combines the 
estimate (jlj) with a generalized version of Schur's Test (Proposition 15.11) . 
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To have a simple idea of the possible applications of Theorems 1, 2 and 3, consider 
the Cauchy problem 

{du 
l m +Hu = 
u(0,x) = Uo(x), 

where H is the Weyl quantization of a quadratic form onK^xR^ (see, e.g., \MW- 
Simple examples are H = — ^A + n\x\ 2 , or H = — ^A — n\x\ 2 (see jl]). The 
solution to ([5]) is a one-parameter family of FIOs: 

u(t, x) = e ltH u , 

with symbol a = 1 and a phase given by a quadratic form $(x, r/), satisfying 
trivially the preceding assumptions ([1]) and (T5]) (see [21] for details). We address 
to Section 7 for a debited study of such operators. 

Finally, Section 6 presents a variant of Theorem 1, cf. fl39l) . and a generaliza- 
tion of Theorem 2, cf. Theorem I6.1[ to the case of FIOs T with symbols in the 
modulation space M 00 ' 1 . This generalizes the known boundedness results on M p of 
pseudodifferential operators with symbols in M 00 ' 1 |25j . and intersects a previous 
result of Boulkhemair [5] on L 2 boundedness of FIOs. We address also to recent 
contribution [8], where the continuity and Schatten-von Neumann properties of 
similar operators when acting on L 2 are proved. 



Notation. We define \t\ 2 = t ■ t, for t G M d , and xy = x ■ y is the scalar product 
on R d . 

The Schwartz class is denoted by <S(M d ), the space of tempered distributions by 
S'(R d ). We use the brackets (/, g) to denote the extension to S(R d ) x S'(R d ) of the 
inner product (/, g) = J f(t)g(t)dt on L 2 (M. d ). The Fourier transform is normalized 
to be f(rj) = Ff{ri) = J f(t)e- 2nitr >dt, the involution g* is g*{t) = g(-t) and the 

inverse Fourier transform is f(r]) = ^ r_1 /(^) = f(—v)- 

Translation and modulation (time and frequency shifts) are defined, respectively, 

by 

T x f(t) = f(t-x) and M v f(t) = e 2 ^ t f(t). 

We have the formulas (T x f)"= M_ x f, {M v fy= T v f, and M V T X = e 2tKixn T x M v . For 
a = («!, . . . , ad), a = (aii, • • • , ctd) G Z^., recall the multi- index notation D a and 
X 13 for the operators of differentiation and multiplication 



d d 
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where t = (ti,...,t d ). We write dx A d£ = ^2j =x dXj A d£j for the canonical 
symplectic 2-form. 

The spaces = l q l^, with weight /x, are the Banach spaces of sequences 
{0"m,n}m,n on some lattice, such that 




(with obvious changes when p = oo or q = oo). 

We denote by Co the space of sequences vanishing at infinity. Throughout the 
paper, we shall use the notation A < B to indicate A < cB for a suitable constant 
c > 0, whereas A x B if c .B < A < cB for a suitable c > 0. 

2. Time-Frequency Methods 

2.1. Short-Time Fourier Transform (STFT). The short-time Fourier trans- 
form (STFT) of a distribution / G «S'(R ) with respect to a non-zero window 
$ G S(R d ) is 

^/(z, 77) = (/, M n T x g) = [ f{t)W^)e- 2mnt dt . 

The STFT V g f is defined on many pairs of Banach spaces. For instance, it maps 
L 2 (R d ) x L 2 (R d ) into L 2 (R 2d ) and x S(R d ) into Furthermore, it can 

be extended to a map from S'(R d ) x S'(R d ) into ^'(M 2 ^). 

We now recall the following inequality [231 Lemma 11.3.3], which is useful when 
one needs to change windows. 

Lemma 2.1. Let g ,gi,l G S(R d ) such that (7,^) ^ and let f G S'(R d ). Then, 
\V g J(x,r])\ < —!—(\V gi f\ * \V gol \)(x, V ), 

for all (x,rj) G R 2d . 

2.2. Modulation Spaces. The modulation space norms are a measure of the 
joint time-frequency distribution of / G S'. For their basic properties we refer, for 
instance, to [231 Ch. 11-13] and the original literature quoted there. 

For the quantitative description of decay properties, we use weight functions 
on the time-frequency plane. In the sequel v will always be a continuous, positive, 
even, submultiplicative weight function (in short, a submultiplicative weight), hence 
v(0) = 1, v(z) = v(—z), and v{z\ + z 2 ) < v(zi)v(z 2 ), for all z,Zi,z 2 G R 2d . 
A positive, weight function \x on R 2d belongs to Ai v , that is, is v-moderate if 
fi(zi + z 2 ) < Cv(zi)fj,(z 2 ) for all z x ,z 2 G R 2d . 
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For our investigation of FIOs we will mostly use the polynomial weights defined 

by 

v s (z) = v s {x, 7]) = (z) s = (1 + \x\ 2 + \ri\ 2 ) s/2 , z = (x, rf) G R 2d . 

Given a non-zero window g G iS(M d ), \i G Ai v , and 1 < p, q < oo, the modulation 
space M p ' q {R d ) consists of all tempered distributions / G <S'(lR d ) such that V g f G 
L p ^ q (M. 2d ) (weighted mixed-norm spaces). The norm on M p,q is 

II/IIm- = \\v 9 f\\ L ™ = ^ iVgffariMwYdxY'dr, 

(with obvious changes when p = oo or q = oo). If p = q, we write M p instead 
of MP' P , and if fi(z) = 1 on R 2d , then we write M p > q and M p for M p > q and M p > p 
respectively. 

Then M p ' q (M. d ) is a Banach space whose definition is independent of the choice 
of the window g. Moreover, if \± G M. v and g G M* \ {0}, then || V^/||x,p.a is an 
equivalent norm for M p ' q (R d ) (see [23, Thm. 11.3.7]): 

II/IIm- x WvjWip. 

2.3. Wiener amalgam spaces. For a detailed treatment we refer to [151 EH1 



Let g G T?(lR 2d ) be a test function that satisfies J2(k i)& 2d T(k,i)9 = 1- Let X(M 2d ) 
be a Banach space of functions invariant under translations and with the property 
that V-X C X, e.g., L p , TU ', or Then the Wiener amalgam space W(X, L^ q ) 
with local component X and global component L^' 9 is defined as the space of all 
functions or distributions for which the norm 

\\f\\w(x,z#*) = ( f {[ \\f ■T( ZliZ2) g\\x) p Kz 1 ,z 2 ) p dz 1 ) q/p dz 2 ^ 
is finite. Equivalently, / G W(X, L^ q ) if and only if 



£(£ll/-%oAMM) p ) 



< OO . 



It can be shown that different choices of g G Z> generate the same space and 
yield equivalent norms. In the sequel we shall use the inclusions relations between 
Wiener amalgam spaces: if B\ B 2 and C\ •—>■ C 2 , 

W(B 1 ,C 1 )^W(B 2 ,C 2 ). 

We now recall the following regularity property of the STFT [9], Lemma 4.1]: 
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Lemma 2.2. Let 1 < p, q < oo, \i G M v . If f G M^ q (R d ) and g G M^,(R d ), then 
V g f G W{J : L 1 ,LPi q ){R 2d ) with norm estimate 

(6) \\Vgf\\ W {rL\L^) < II / II Ik II A/i • 

We also give a slight generalization of [23, Proposition 11.1.4] and its subsequent 
Remark. 

Proposition 2.1. Let X be a separated sampling set in R 2d , that is, there exists 
5 > 0, such that ia£ x ,yeX:x^y \ x ~ y\ > Then there exists a constant C > 
such that, if F G W(L°°,L^ q ) is any function everywhere defined on R 2d and lower 
semi- continuous, then the restriction F\ x is in £ p ~ 9 , where jl = fi\ x , and 

\\F\ X \\^ 9 < C\\F\\ W ( LO oj$<iy 

Proof. One uses the arguments of [23], Proposition 11.1.4] and its subsequent Re- 
mark. We just shall highlight the key points that make those arguments to work 
under our assumptions. 

First of all, if (r, s) G Z M and X is separated, then the number of sampling 
points of X in (r, s) + [0, l] 2d is bounded independently of (r, s). 

Secondly, forx&X such that x G (r,s) + [0, l] 2d , 

\F(x)\fi(x) < C?||F-T( r> ,)X[ ,i]a«t|U»//(r,s), 

Indeed, since F is everywhere defined on R 2d and lower semi-continuous we have 
SU P \ f\ — esssup|/| on every box, whereas fi(x) < C/x(r, s) because \i is w-moderate 
and v is bounded on [0, l) 2d (cf. (23J Lemma 11.1.1]). □ 

2.4. Gabor frames. Fix a function g G L 2 (R d ) and a lattice A = aL d x for 
a, (3 > 0. For (m, n) G A, define g mjn := M n T m g. The set of time-frequency shifts 
C/(<7, a, /3) = {g m>n , (m, n) G A} is called Gabor system. Associated to Q(g, a, (5) we 
define the coefficient operator C g , which maps functions to sequences as follows: 

(7) (C g f) m ,n = (Cff) m ,n ■= (/, 9m,n), (m, 7l) G A, 

the synthesis operator 

-Dp c ^ ^ c m ^ n T m M n g , c {'-"m,n}(m,n)EA 

(m,n)GA 

and the Gabor frame operator 

(8) S 9 f = Sf f := D B S B f = (f,9m»)9m,n- 

(m,n)eA 

The set Q(g,a,/3) is called a Gabor frame for the Hilbert space L 2 (R d ) if 5^ is 
a bounded and invertible operator on L 2 (R d ). Equivalently, C g is bounded from 
L 2 (R d ) to Z 2 (aZ d x pZ d ) with closed range, i.e., ||/|| L 2 x \\C g f\\p. If G(g,a,{3) is a 
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Gabor frame for L 2 (R d ), then the so-called dual window 7 = S~ l g is well-defined 
and the set a, 0) is a frame (the so-called canonical dual frame of Q(g, a, (3)). 
Every / G L 2 (M. d ) posseses the frame expansion 

(9) / = (f'9m,n)lm,n= (f ' ~fm,n) 9m,n 

(m,n)6A (m,n)£A 

with unconditional convergence in L 2 (M. d ), and norm equivalence: 

\\f\\ L , x \\C g f\\ P x ||C 7 /|| J2 . 

This result is contained in [23, Proposition 5.2.1]. In particular, if 7 = g and 
1 1 g 1 1 £2 = 1 the frame is called normalized tight Gabor frame and the expansion (Q 
reduces to 

(10) /= ^ (f>9m,n)g m ,n- 

(m,n)£A 

If we ask for more regularity on the window g, then the previous result can be 
extended to suitable Banach spaces, as shown below [191 26J. 

Theorem 2.2. Let \i G M. v , Q(g,a,/3) be a normalized tight Gabor frame for 
L 2 (R d ), with lattice A = aZ d x (3Z d , and geM$. Define A = AV 
(i) For every 1 < p, q < 00, C g : M^' q —>■ l p ~ q and D g : — > countinuously 
and, if f G M^ ,q , then the Gabor expansions ({TO]) converge unconditionally in M^ ,q 
for 1 < p, q < 00 and all weight fi, and weak*-M^° unconditionally if p = 00 or 
q = 00. 

fiij T/ie following norms are equivalent on M v ^ q : 

(11) 11/11*™ - l|c s /lb r . 

We also establish the following properties. Denote by the closure of the 

Schwartz class in M^ ,q . Hence, M p,q = M^' q if p < 00 and g < 00. Also, denote by 

the closure of the space of eventually zero sequences in Z~' 9 . Hence = if 
p < 00 and q < 00. 

Theorem 2.3. Under the assumptions of Theorem \2.S\ r for every 1 < p, q < 00 the 
operator C g is continuous from M^ ,q into Z?' 9 , whereas the operator D g is continuous 
from l p ~' q into Mf; q . 

Proof. Since C g is continuous from M v ^ q into I™ it suffices to verify that, if / is a 

Schwartz function then C g (f) G l^ 9 . This follows from the fact that C g (f) G Z~. 
Similarly, for D g it suffices to verify that, if c is any eventually zero sequence, then 
D g(c) G MP' q . This is true because D g (c) G M*. □ 
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3. Almost diagonalization of FIOs 

For a given function / on R d the FIO T with symbol a and phase $ can be 
formally defined by 

(12) Tf(x)= [ e 2 ^^a(x, V )f( V )d V . 



To avoid technicalities we take / G iS(lR d ) or, more generally, / G M 1 . If a G L°° 
and the phase $ is real, the integral converges absolutely and defines a function in 
L°°. 

Assume that the phase function 77) fulfills the following properties: 

(i) $ G C°°(R M ); 

(ii) for z = (x, 77), 

(13) \d2$(z)\<C a , \a\>2; 

(iii) there exists 5 > such that 

(14) I det d 2 x ^(x,r])\ > 5. 
If we set 

y = V v ®(x,rj) 



(15) 



f = V x $(z,»7), 



and solve with respect to (x, £), we obtain a mapping x, defined by (x,£) = x{Vi v)y 
which is a smooth bilipschitz canonical transformation. This means that 

- x is a smooth diffeomorphism on R M ; 

- both x and x~ l are Lipschitz continuous; 

- x preserves the symplectic form, i.e., 

dx A = dy A g?77. 

Indeed, under the above assumptions, the global inversion function theorem (see 
e.g. [32J) allows us to solve the first equation in (fT5l) with respect to x, and 
substituting in the second equation yields the smooth map x- The bounds on the 
derivatives of x, which give the Lipschitz continuity, follow from the expression for 
the derivatives of an inverse function combined with the bounds in (ii) and (iii). 
The symplectic nature of the map x is classical, see e.g. [7]. Similarly, solving the 
second equation in (fT5|) with respect to 77 one obtains the function x _1 with the 
desired properties. 

In this section we prove an almost diagonalization result for FIOs as above, with 
respect to a Gabor frame. Here we consider the case of regular symbols. In Section 
|6]we will study the case of symbols in modulation spaces. 
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Precisely, for a given N e N, we consider symbols a on M? d satisfying, for 

z = (x,ri), 

(16) \dy(z)\ < C a , a.e. C G R 2d , \a\ < 2N, 

here <9" denotes distributional derivatives. 

Our goal is to study the decay properties of the matrix of the FIO T with respect 
to a Gabor frame. For simplicity, we consider a normalized tight frame G(g,a,/3), 
with g e S(R d ). 

Theorem 3.1. Consider a phase function satisfying (i) and (ii) and a symbol 
satisfying f flBT) . There exists Cn > such that 

(17) \{Tg m ,n, 9m',n'}\ < Cjv(V*$(m', n) - (n',m))- 2N . 

Proof. Recall that the time-frequency shifts interchange under the action of the 
Fourier transform: (T x f) A = M_ x f and (M rj f) A = T v f, besides they fulfill the 
commutation relations T X M V = e~ 2mxr> M V T X . Using this properties, we can write 



Tg m ^{x)M n ,T m ,g(x) dx 

e 2m ^a(x } rj)T n M- m g{rj)M- n ,T m ,g{x) dxd V 
M(p,- m )T(p- n ) (e^^aix^)) g( V )M- n ,T m ,g(x) dxd v 
T(- m ^)M_ (< o)M (0) _ m) T (0 ,_ n ) (e 2m ^a(x, rj)) g(x)g(rj) dxd V 

e 2 m Mx+m', V+ n)-(n',m)ix+ m ', V )] a ^ x + m ' ^ + n)g(x)g(rj) dxdf] 



Since $ is smooth, we expand $>(x,r]) into a Taylor series around (m',n) and 
obtain 

$(x + m! , T] + n) = $(m', n) + V 2 $(m', n) ■ (x, rf) + $2,(m',n)( x i v) 
where the remainder is given by 

(18) ®2,(m',n)(x,v) = 2 V / (1 - t)<9 Q $((m» +t(x, V )) dt^^. 

Whence, we can write 
\{Tg m,ni 9m', n' / \ 

e 2,rf{[V.*(m',n)-(n',m)].( ! r,,)} e 2 9 r<* 3i(m , in) (^) (T ( a . + m > ^ + n)g{x)g{rf) dxdf] 
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For iVeN, using the identity: 

(l-A^e 2 ™^*^^"''"*^)} = (27T(V,$(m', n)-(ri, m )))2^ e 2^{[V,*(m',n)-(n',„ 

we integrate by parts and obtain 
\{Tg m,m 9m', n' / \ 

1 



2-Ki{[V z $(m',n)-(n',m)]-(x,T))} 



(27r(V 2 <£>(m',n) - {n',m))) 2N 

x (1 - A/ \ e 2 ™*2,<, m >, n )(x,v) a ( x + m ' }V + n)g(x)g{r])] dxdr] 

By means of Leibniz's formula the factor 

(1 - A/ \e 2wi ^'^ {x ' v) a{x + m', V + n)g(x)g{r)) 
can be expressed as 

|a|+]^|+|7l<2JV 

where p(d^Q 2 ,(m' ,n))( z ) is a polynomial made of derivatives of $2,(m',n) of order at 
most \a\. 

As a consequence of (ii) we have d"^2,(m',n)(z) = 0((z) 2 ), which combined with 
the assumption ffTB]) and the hypothesis g G 5(R d ) yields the desired estimate. □ 

Remark 3.2. More generally one can consider symbols satisfying estimates of the 
form 

(19) \dy(z)\ < C a fi(z), a.e. z G R 2d , \a\ < 2N, 

with n G A4 V and also more general windows g. Indeed, by arguing as above and 
using 

\d^a(z+(m',n))\ < C^(z + (m',n)) < C'C p v(z)ii(m' , n), 
one deduces the decay estimates 

/i(m', n) 



m,n> Qm' ,n' ) \ < Cjv 



(V 2 $(m', n) — (V, m)) 



2V : 



provided the integral 



(20) 



QA 7 Md H $2,( ? <™))(^M^)^ ® g){i 

|o| + l/3| + l7l<2V 



converges. This is guaranteed if, e.g., (z) 2N v(z)d](g ® g)(z) G L x . 
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We now assume the additional hypothesis (iii) on the phase, and rewrite ( jTTl) in 
a form convenient for the applications to the continuity of FIOs in the next section. 
We need the following lemma. 

Lemma 3.1. Consider a phase function $ satisfying (i) ; (ii), and (iii). Then 

(21) | V x $(m', n) — n'\ + | V^$(m', n) — m\ > \x(m, n) — m'\ + \C,{fn, n) — n'\, 
where (y,i]) i — > (x, £) is t/ie canonical transformation generated by $. 

Proof. It suffices to prove the following inequalities: 

(22) | V v &(m', n) — m\ > |a;(m, n) —m'\, 

(23) (V^r^n) -n'| > |£(m,n) -n'| - C| V r ,$(m / , n) -m|. 
We observe that, by (|15l) . we have 

(24) J/ = V„$(a;(y^)^) V(y,7?)eR 2d 
and 

(25) V x $(x,7?) =f(V„$(x, 77), 77) V( I)t) )6M M . 
Hence, we have m = W v (x(m, n), n), so that 

|V^$(m', n) — m\ — \V v §(m' , n) — Vr,$(x(m, n), n)| 
> \x(m, n) — m'|, 

wehere the last inequality follows from the fact that, for every fixed 77, the map 
x 1 — > VQ v {x,v) has a Lipschitz inverse, with Lipschitz constant uniform with 
respect to r\. 

This proves (1221) . 

In order to prove (T231 we observe that, in view of (1251) . it turns out 

V a; $(m / , n) -n' = ^(V^$(m', n), n) - 

= £(m + V^$(m', n) — m,n) — n 

= C,(m, n) — n' + 0(V^$(m', n) — m)). 

where the last inequality follows from the Taylor formula for the function y 1 — > 
£,(y,n), taking into account that the function £ has bounded derivatives. 

This proves ff23|) . □ 

Combining the previous lemma with ffTTl) we obtain the following result. 

Theorem 3.3. Consider a phase function $ satisfying (i), (ii), and (iii), and a 
symbol satisfying f fTBT) . Let g G 5(lR d ). There exists a constant > such that 

(26) \(Tg m ,n,g m >,n>)\ < C N ( X (m,n) - (m',n'))~ 2N , 
where x is the canonical transformation generated by $. 
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This result shows that the matrix representation of a FIO with respect a Gabor 
frame is well-organized, similarly to the results recently obtained by 0, [27J in 
terms of shearlets and curvelets frames. More precisely, if a G S® Q , namely if (fTBl) 
is satisfied for every N G N, then the Gabor matrix of T is highly concentrated 
along the graph of x- 

4. Continuity of FIOs on Ml 

In this section we study the continuity of FIOs on the modulation spaces as- 
sociated with a weight function fx G -M„ 3 , s > 0. We need the following preliminary 
lemma. 

Lemma 4.1. Consider a lattice A and an operator K defined on sequences as 

(Kc) x = 2J K x ,yC v , 

where 

sup 2_] \K\A < °°; SU P / J \K\A < °°- 

ASA ygA 

TTien ii" zs continuous on / P (A) /or even/ 1 < p < oo and moreover maps the space 
Co(A) of sequences vanishing at infinity into itself. 

Proof. The first part is the classical Schur's test (see e.g. [23| Lemma 6.2.1]). The 
second part follows in this way. Since we know that K is continuous on l°° and the 
space of eventually zero sequences is dense in Co, it suffices to verify that K maps 
every eventually zero sequence in cq. This follows from the fact that any eventually 
zero sequence belongs to I 1 and therefore, since K is continuous on I 1 , is mapped 
in l l <— > Co. □ 

We can now state our result. 

Theorem 4.1. Consider a phase function satisfying (i) , (ii), and (iii), and a symbol 
satisfying ({16]) . Let < s < 2N — 2d, and ji G M. Vs - For every 1 < p < oo, T 
extends to a continuous operator from M? ox into M?, and for p = oo it extends to 
a continuous operator from M^ x into . 

Recall that is the closure of 5(IR d ) in M^°. Moreover, observe that /i o % g 
Ai Vg . Indeed, v s o \ x v s , due to the bilipschitz property of x- 

Proof. We first prove that 

l|77|| M * ox < CH/IIm., 

for every / G iS(IR d ). This proves the theorem in the case p < oo, since <S(M d ) is 
dense in Mf.. 
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We see at once that, since a G L°°, T defines a bounded operator from M 1 into 
Lo° ^ M oo_ Hence, for all / G S{R d ), we have Tf G M°° and Theorem O shows 
that ||/|| m p ox > \\C g (f)\\ lUx and \\Tf\\ M » < \\C a {Tf)\\%- ° n the other hand, the 
expansion fllOp holds for / with convergence in M 1 . Therefore 



Tf = ^2(f,g m ,n)Tg ri 



with convergence in M°°. Hence, 
Cg{Tf) m ',n' = (Tf,g m/>n ,) = ^2(Tg 

m,ni 9m', n' )(f,9m,n) = ^2( T 9 m,ni 9m', n' )C g (f) m,n- 
m,n m,n 

Therefore we are reduced to proving that the matrix operator 

(27) {c m rl J' I > ^ (7 9m,ni 9m' '^m,n 

is bounded from — into Z~. This follows from Schur's test (Xemma 14.11) if we 
prove that, upon setting 

^■m' ,n' ,m,n (Tg m,ni 9m', n' 

we have 

(28) K-m' ,n' ,m,n G l m ,rj'm',n'i 

and 

(29) Km' ,n' ,m,n G l m ' ,n'^m,n- 

In view of fl26l) we have 

(30) |# mW , m , n | < (x(m, n) - (m', n')>~ 2iV+s ^ , f^'^ ; , ^ 

(X(m, n) - (m', n')) s p(x{m, n)) 

Now, the last quotient in ( [301) is bounded because \i is t> s -moderate, so we deduce 

(EHD- 

Finally, since x is a bilipschitz function we have 

(31) - (m',n')| x |(m,n) - ^(m'X)! 

so that f)29p follows as well. 

The case p = oo follows analogously by using Theorem 12.31 (with p = q = oo), 
and the last part of the statement of Lemma 14. 1 1 □ 

Remark 4.2. Theorem 14.11 with v = 1 gives, in particular, continuity on the 
unweighted modulation spaces M p . If moreover p = 2, we recapture the classical 
L 2 -continuity result by Asada and Fujiwara pQ. 

Also, Theorem 14. II applies to p = v t , with \t\ < s. In that case we obtain continuity 
on M Vt , because v t ° x x v t- 
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5. Continuity of FIOs on M p ' q 

In this section we study the continuity of FIOs on modulation spaces M p,q possi- 
bly with p q. As shown in Section under the assumptions of Theorem 14. II such 
operators may fail to be bounded when p ^ q. The counterexample is given by the 
phase $(or, rf) = xrj + \x\ 2 /2, and symbol a — 1, which does not yield a bounded 
operator on M p ' q , except for the case p = q. Here the obstruction is essentially due 
to the fact that the map x i — > S7 x <&(x,r)) has unbounded range. Indeed we will 
show, for general phases, that if such a map has range of finite diameter, uniformly 
with respect to 77, then the corresponding operator is bounded on all M p,q . To this 
end we need the following result. 

Proposition 5.1. Consider an operator defined on sequences on the lattice A = 
aZ d x pZ d by 

(^Kcj m '^ n i ^ K- m i n i m n c m n . 



(1) If K e ln l n' l m' l ln> K is continuous on l\V£. 
(ii) If K e l^ln^m'} K ^ s continuous on l^l^. 

(Hi) IfKe Wlrl^t n and moreover K G C',X,n H C.X'.n" the 

operator K is continuous on l p,q = l^l^ for every 1 < p, q < 00. 

(iv) Assume the hypotheses in (Hi). Then K is continuous on all l p,q , 1 < p, q < 00. 

Recall that l p,q is the closure of the space of eventually zero sequences in l p,q . 

Proof, (i) We have 

||-K"c||,i ,oo < Vsup V"|X„ 



L m' ,rv ,m,n 



, m' 
n' m,n 



\ 1 m' 

n \ n m 

ri n i ™/ 111 



SUp \C m ,r 



(ii) It turns out 



,m,n \'-"m,n 



<supV V |-KTm>' 

m/ \m,n 

< SUp SUp ^ |i^m',n',m,n| ^ ' 

n' \ m , / 

n \ m! m / 



< ll^llz-te*i INI*-/ 1 - 



(m) Since the statement holds for p = q by the classical Schur's test, and for 
(p, q) = (1, 00) and (p, q) = (00, 1) by the items (2) and (ii), it follows by complex 
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interpolation (see (3) on page 128 and (15) on page 134 of [37]) that it holds for 
all (p, q), except possibly in the cases q = oo, 1 < p < oo. For these cases we argue 
by duality as follows. 

In order to prove the continuity of K on l°°l p , it suffices to verify that for any 
sequences c = (c mj „) G and d = (d m ',n') G J^m'' with d m / jn / > 0, we have 

(32) ^ \( Kc )m',n'\d m >,n' < || c ||j<?#, 

m' ,n' 

Now 

^ ^ | (-^^)m',n' I d"m r ,n' — ^ ^ ^ ^ \^m' ,n' ,m,n\\^m,n\dm' ,n' 
m',n' m',n' m,n 

^ ^ J ^ ^ |-^m' ,n' ,m,n\dm' ,n' J |Cm,ri| 
m,n \m',n' / 

where is the operator with matrix kernel K m ,n,m',n' = \K m ',n',m,n\- Since it 
satisfies the same assumptions as K, it is continuous on l 1 !? , which gives (1521 . 
(if) Since X is continuous on l p,q and by the definition of l p,q , it suffices to verify 
that K maps every eventually zero sequence in l p,q . This follows from the fact that 
K maps every eventually zero sequence in I 1 Z 13 ' 9 , because K is bounded on Z 1 . □ 

Theorem 5.2. Consider a phase function $ satisfying (i), (ii), and (iii), and a 
symbol satisfying ( 1161) . u^/j N > d. Suppose, in addition, that 

(33) sup |V x $(a;,77) - V x $(:r', ??)| < oo. 

x,x' ,r) 

Then the corresponding Fourier integral operator T extends to a bounded operator 
on M p,q for every 1 < p, q < oo and on M p,q if p = oo or q = oo. 

Proof. By arguing as in the proof of Theorem 14.11 it suffices to prove the continuity 
on l p ' q = l^lfn if p < oo and q < oo, or l p,q if p = oo or q = oo, of the operator 

{c m ,„} — £ T w> , m>n c Win; 

where 

T m ',n',m,n \T ' 9m,ni 9m' ,n') • 

By applying 15.11 it suffices to verify that 

(34) {Tm'^'^m^} G Z n l n ,l m il m , 

(35) {^m'^m^} G l n il n l m l m , , 
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because we already see from ([26]) and (JSTJ) that {T m ^ n ^ myn } e lm,Jm',n> n ^m',n'^m,n- 
Let us now prove AM]) . It follows from (|T7|) and (123"j) that 

|T m / in ^ m , n | < ( 1 + \V x $(m , n) - n I + |V,,$(m , n) - m| 
< ^1 + |£(m, n) — r/| 2 + | V^$(m', n) — m 

By ( TT5l) we have 



2 , 



< (1 + |£(m, n) - n'\) N (1 + | V„$(m', n) - m|) _iV . 



£(</, 77) = V a $(x(y, 77), 77), V{y, 77) G M M , 
so that the hypothesis (1331) yields 

£(m,n)=£(0,n) + O(l). 

Hence ( 1341) follows. 

We now prove As above, it follows from f[26"j) and (1221) that 

l^m'y.m.nl ~ ( x + |V a .$(m', n) - n'| 2 + |V„$(m', n) - m| 2 ) 

-A 7 " 



< ^1 + |V x $(m', n) — n'| 2 + |a;(m, n) — m'| 2 j 

< (1 + |V x $(m', ra) - n'])"* (1 + |z(m, n) - m'|) _JV 

By ( ]33l we have 

V x $(m',n) = V*$(O,n)+0(l), 

so that 

1 + I V x $(m', n) - n'\ > 1 + | V x $(0, n) - n'| 
(36) >l + \n-^(n% 

where ip is the inverse function of the bilipschitz function 77 i — ► V x $(0, rj). There- 



fore we obtain ( l35lh 

This concludes the proof. □ 

Example 5.3. Theorem I5.2I applies, in particular, to phases of the type 

rj) = xrj + a(x, rj) where \d%d!?a(x,r})\ < C a ,p for every 2\a\ + > 2. 

In the special case when a(x,rj) = a(rj) is independent of x and the symbol a = 1, 
the FIO reduces to a Fourier multiplier 

Tf(x)= [ e^ ixr >e 27ria ^f(r])dr] 



and we reobtain the result of [2j Theorem 5] on the continuity of T on all M p,q , 
1 < P, g < 00. 
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6. Modulation spaces as symbol classes 

In what follows we shall rephrase the quantity \(^T g m ,ni 9m> ,n')\ in terms of the 
STFT of the symbol a, without assuming the existence of derivatives of a. This will 
be applied to prove the continuity of FIOs with symbols in M 00 ' 1 on modulation 
spaces M p . 

The same arguments as in Theorem 13.11 yield the equality 

(Tg m , n , g m >,n>) = e 2 ™ / e 2m ^a(x, r ? )M ( _ ri ,,_ m) T (m ,, n) (^ ® g)(x, rj) dxd V . 
Expanding the phase $ into a Taylor series around (m',n) we obtain 

rj) = $(m', n) + V z $(m', n) ■ (x - m', rj - n) + T {m ^ n) $ 2 ,{ m ' ,n){x, v) 

where the remainder $2,(m',n) is given by ([TBI . 

Inserting this expansion in the integrals above, we can write 

(37) (Tg m>n ,g m > >nl ) = e ^(»nn+*("» , .»)-V,*(m',n).(m',n)) f ^iV^m' ,n)( X , V ) 

JR d 

x M ( _ nVm) T (m; , n) e 2 ^, (m '.„)(^)^ g, ^ ^ 

Defining 

(38) *7) := e 2 ^^'.")^^®^)^^), 

and computing the modulus of the left-hand side of (1371) . we are led to 

(39) \(Tgm, n ,g m '<n') \ = |V* (m , n) e((m', n), (n - V x $(m', n), m - V„$(m', n)))| 

Observe that the window ^(m^n) of the STFT above depends on the pair (m', n). 
We now study the continuity problem of T when the symbol a is in the modu- 
lation space M 00 ' 1 . 

Theorem 6.1. Consider a phase function satisfying (i) , (ii), and (in), and a symbol 
o G M 00,1 . For even/ 1 < p < oo, T extends to a continuous operator on M p , and 
for p = oo it extends to a continuous operator on M°° . 

By arguing as in the proof of Theorem 14.1^ it suffices to prove the continuity on 
V if 1 < p < oo and on l°° = cq, of the operator (|27p . In wiew of Schur's test 
(Lemma 14. ip and (1391) . it suffices to prove the following result. 

Proposition 6.2. Consider a phase function $ satisfying (i) and (ii) and (iii) and 

a symbol a G M 00,1 . If we set 
(40) 

Zm,n,m',ri '■= {{ m ' i n )i ( n '~ Va;$(m', n), m— V^$(m', n))), m,m! G aZ d , n, n' G /3Z d , 
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then, 

(41) SUp \ V * { m>. 
(m,n)eA {mly)eA 

(42) sup \ V * (m '. 

( m '"'' 6A Kn)eA 

We need the following lemma. 

Lemma 6.1. Let * G S(IR M ) raf/i ||\I/o||l 2 = 1 and \I/( m » be defined by fl38l) . 
wt/i (m',n) G A = aZ d x /3Z rf , and # G «S(R d ). T/jen, 

(43) / sup \Vv Jtf (w)\dw<oo. 

JK 4d (m',n)eA 

Proof of Lemma \6.1\ We shall show that 

(44) l^,n)*oH| < C(w)-^ d+1 \ V(m» G A. 
Using the switching property of the STFT: 



( Vf g)(x,r ] ) = e- 2 ^(V g f)(-x,-r ] ), 

we observe that \V^ {m , ^ (w 1 ,w 2 )\ = \V 9o ^ ^ ^(-Wx, -w 2 ) , and by the even 
property of of the weight (•), relation (144]) is equivalent to 

(45) \V^ (m r, n) {w)\ < C(w)- {M+l \ V(m',n) G A. 

Now, the mapping V^ is continuous from S{R 2d ) to 5(M 4d ) (see [231 Chap. 11]). 
This means that there exists M G N, K > 0, such that 

(w) M+1 



< 1 1 Vy^tm' ,n){-) |U°°(R«)(w) 
<A^ ||^ 7 ^(m',n)||L-(K^)(«;) 



-(4d+l) 



| 7 |+|5|<Af 



for every (m',n) G A. We now claim that ^/( m ',n) G 5(1R M ) uniformly with respect 
to (m', n). This is proved as follows: the function e 27n * 2 -( m ' j s j n C°°(R 2d ) an d 
possesses derivatives dominated by powers ((x,r])) k , fceN, uniformly with respect 
to (m',n), due to flE]); since (5 <g> <?) G <S(M M ), it follows that ^ (m » G 5(M 2d ), 
with semi-norms uniformly bounded: 

PM(*(m',n)) := ^ ||# 7 -X'**(m',n)|U<» < C M , V(m', 7l) G A. 

|7|+|5|<M 



20 ELENA CORDERO, FABIO NICOLA AND LUIGI RODINO 

Consequently, 

\V^ {m/>n) (w)\ < K J2 P 7 ^^( m ',n)IU-(R-)(^)- (4d+1) < KC M (w)-( 4d+1 \ 

\j\+\S\<M 

for every (m', n) G A, as desired. □ 

Proof of Proposition WB . We shall prove (HI]) . First, Lemma [2. II for gi = 7 = 
yields 

\V* (m , tn) <r(z)\ < {\V 9o a\ * |VV n) * |)(*), z G R 4d , 

so that 



I Zm.n.m' ,n' 



^(m'^i 

/ ^2 \ V ^0 a ( Z m,n,m',n' ~ Vj)\ | Vq> (m , n) ^q(w) \ dw 

^ R4d (m',n')eA 

/ Y] \Vy (r{z m>n>m > >n i - w)\ sup \V 9 ^ (w)\dw 

JR 4d , 7A cA (m',n)eA 

< SUp \Vvo^(Xm,n,m',n' / SUp \V 9 , W (w)\dw 

wm 4d / i , w JM. 4d (m',n)eA 



(m',n')eA 

< 
< 



(m',n')eA 

< C SUp ^ \V*o a { z m,n,m> >' 
[m ,n )£A 

where the last majorization is due to Lemma [6.11 Since, 

Zm,n,m' ,n' 

— w)\< 2j SU P |V* o a(Mi,Z m 

,n,m' ,n' ,W2 ) \ 5 

(m',n')eA (m',n')eA uieRM 

with 

Zm,n,m» 2 : = ( n ' - V x $(m', n),m- V ri $(m', n)) - w 2 , w = (w u w 2 ) G 
we shall prove that 



(46) E SU P |V^ cr(tii, fm, n ,m» 2 ) | ^ IM 



( m ',nOeA«i eR2d 



uniformly with respect to (m, n) G A, w 2 G M M For every fixed (m, n), the set 

# = ^TO,n,i«2> g iven b y 

Xm,n,w 2 — {Zm,n,m' ,n' ,w 2 ', {m',n') G A, W 2 G R M }, 

is separated, uniformly with respect to (m, n), u>2- Indeed, given (m^ni) 7^ 
(m' 2 ,n' 2 ), if mi 7^ m' 2 , 

l^n.n.mi.ni.toa ~ 2m,n im » 2 | > |V„$(mi,7l) - V„$(m£,7l)| > C|mi - m' 2 | > aC, 
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uniformly with respect to (m, n), io 2 , because the mapping x i — > V, 7 $(x, 77) has an 
inverse that is Lipschitz continuous, thanks to ([TBI and (|14p . On the other hand, 
if m[ = m 2i 

Hence, <Y is separated uniformly with respect to (m, n), W2- Now, we apply Propo- 
sition 12.11 (with p — q = 1) to the function 

F(u 2 ) = sup |V* cr(ui, w 2 )|, «2 £ K M , 

which is lower semi-continuous, being Vq, a continuous. We obtain 

(47) II^Ull^ 1 < Cll^Hty^.i 1 ) = II Vt>o a \\w(c,L°°' 1 )- 

If the symbol cr is in M 00 ' 1 , by Lemma O the STFT Vy a belongs to the Wiener 
amalgam space W{TL X , L 00 ' 1 ), and 

ll^o^llw/^L 00 ' 1 ) ^ II^o^IIh/^L 1 ,!, 00 ' 1 ) ^ HcIIa/ 00 ' 1 II^oIIm 1 - 

The first inequality is due to TL X <^-> C and the inclusion relations between Wiener 
amalgam spaces. Combining this inequality with (|47|) we obtain (H6l) . uniformly 
with respect to (m,n) and W2, that is (HIT) . 

The estimate (j4"21) is obtained by similar arguments. □ 

Remark 6.3. We observe that the continuity on M 2 = L 2 of FIOs as above, with 
symbols in M 00 ' 1 , was already proved in [5] by other methods. 

7. The case of quadratic phases: metaplectic operators 

In this section we briefly discuss the particular case of quadratic phases, namely 
phases of the type 

1 1 

(48) $(x, n) = -Ax ■ x + Bx ■ 77 + ' V + Vo ■ x — xo ■ r], 

where Xo, t]q G M. d , A, C are real symmetric d x d matrices and B is a real d x d 
nondegenerate matrix. 

It is easy to see that, if we take the symbol a = 1 and the phase PHI) , the 
corresponding FIO T is (up to a constant factor) a metaplectic operator. This can 
be seen by means of the easily verified factorization 

(49) T = M^V A D B T- X V C TT X ^ 

where U a and Uc are the multiplication operators by e mAx ' x and e 7nCr? ' ?? respectively, 
and Db is the dilation operator / 1— > f(B-). Each of the factors is (up to a constant 
factor) a metaplectic operator (see e.g. the proof of [311 Theorem 18.5.9]), so T is. 
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The corresponding canonical map, defined by ( TToT) . is now an affine symplectic 
map. For the benefit of the reader, some important special cases are detailed in 
the table below. 



operator 


phase rf) 


canonical transformation 


T 


(x -Xo)-T) 


x(y,v) = (y + xo,v) 


M m 


(V + Tjo) ■ X 


x(y,v) = (y,v + vo) 


D B 


Bx ■ 7] 




U a 


X ■ T) + \Ax ■ X 


x(y,v) = (y,v + Ax ) 



However one should observe that there are metaplectic operators, as the Fourier 
transform, which cannot be expressed as FIOs of the type (fT2|) . 

Metaplectic operators are known to be bounded on M£ , see e.g. [23] Proposition 
12.1.3]. This also follows from Theorem l4.ll Indeed, since x is a bilipschitz function, 
we have v s ° X x v s- 

Also, Theorem 15.21 applies to quadratic phases whose affine symplectic map x 
is (up to translations on the phase space) defined by an upper-triangular matrix, 
which happens precisely when A = 0. Indeed, we obtain the map x by solving 

y = Bx + Cr] — x 

£ = Ax + Br] + r] . 

The phase condition (fill) here becomes 

det B ^ 0, 



so that B is an invertible matrix and x 
mapping x ■ (y, v) 1 — ► (x, is given by 

B- 1 -B- l C 



x 




a. 





AB- 



B — AB~ X C 



B~ x y - B~ x Cri + B~ x x Q . Whence, the 



B- l x Q 





y 


+ 









AB^xo + rj 



When A = the phase $ satisfies f[33|) and, consequently, the corresponding opera- 
tors are bounded on all M p,q . This can also be verified by means of the factorization 
(jl9|) (with A = 0). Indeed the continuity of the operators M Vo , T Xo and Db is eas- 
ily seen, whereas that of the Fourier multiplier T^UqF was shown, e.g., in [25| 
Lemma 2.1]. 

On the other hand, generally the metaplectic operators are not bounded on M p,q 
if V 7^ Q- An example is given by the Fourier transform itself (see [2]). An example 
which instead falls in the class of FIOs considered here is the following one. 



Proposition 7.1. The multiplication Uj d 
p,q < oo, with p 7^ q. 



is unbounded on M p,g , for every 1 < 
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Proof. We have U[ d f(x) = e m ' x ' f(x). For A > 0, we consider the one-parameter 
family of functions f(x) = e ~ nX W 2 e M p ' q , so that f(r)) = \-d/-i e -n(i/\)\ v \\ For 
every 1 < p, q < oo, by [TH Lemma 5.3], we have 

(A + l) rf( H) 

Since Uf^e'^-^ 2 , so that Uf(r]) = (A - i)-^ e -^V(^i)\v\\ the same formula 
as above yields 

[( A + 1)2 + i]KH) 



11^/11 



A^(A 2 + A + 1)2 ( F-^) 
As A — > 0, we have 

so that, if we assume ||?7/||mp.9 < then 1/p — 1/q > 0, that is p < q. 

Moreover, the same argument applies to the adjoint operator U*f(x) = e~ m ^ f(x). 

Now we show that p = q. By contradiction, if U were bounded on M p,q , with 
p < q, its adjoint U* would satisfy 

\\U*f\\ MP >, q ><C\\f\\ MP ,, q ,, V/G<S(R d ), 

with q' < p', which is a contradiction to what just proved. □ 
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